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Abstract 

The detectability of planetesimal impacts on imaged exoplanets can be measured using Jupiter during 
the 1994 comet Shoemaker-Levy 9 events as a proxy. By integrating the whole planet flux with and without 
impact spots, the effect of the impacts at wavelengths from 2-4 ftm is revealed. Jupiter’s reflected light 
spectrum in the near-infrared is dominated by its methane opacity including a deep band at 2.3 /rm. After 
the impact, sunlight that would have normally been absorbed by the large amount of methane in Jupiter’s 
atmosphere was instead reflected by the cometary material from the impacts. As a result, at 2.3 /im, where 
the planet would normally have low reflectivity, it brightened substantially and stayed brighter for at least 
a month. 
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1. Introduction 

The frequencies with which giant planets and cold circumstellar debris disks encircle old (i.e. >1 Gyr) 
Sun-like stars are both ^20% (??). The debris disks are generated in the collisions and evaporation of 
planetesimals, analogous to the comets and asteroids of today’s Solar System. That most debris disks 
contain cold, ^^50 K, dust indicates that the reservoirs of parent bodies are at '^100 AU from the parent 
stars, in a region analogous to our Edgeworth-Kuiper Belt. These cold debris disks appear to be equally 
common for planet-hosting and non-planet hosting stars (?), though the known planets in these systems are 
all within a AU of their stars. Detectable disks are at least 10-100x dustier than our Edgeworth-Kuiper 
Belt (??), with correspondingly larger parent body populations and are thus even more subject to frequent 
collisions between bodies. 

In today’s Solar System, planetesimals infrequently hit planets; although, in the past, such impacts were 
much more common. The cratering rates on the Moon and Mars indicate that the impact rates were high 
prior to ~800 Myr after the formation of the Solar System, perhaps ending in a “Late Heavy Bombardment” 
(see e.g. ?, and references therein). Today, a comet larger than 1 km in diameter impacts Jupiter only every 
few hundred years (???). One such collision was that of Jupiter and Shoemaker-Levy 9 (SL9) in July 1994. 
SL9 was a Jupiter Family comet, a class of objects believed to originate in the Edgeworth-Kuiper Belt (?). 
During a pass of Jupiter in 1992 the comet was tidally disrupted and broke into ^25 fragments (?). Three 
of these fragments were still quite large - >1 km in diameter (??). Fifteen fragments produced impacts 
detected by the Hubble Space Telescope (?), as the comet impacted the planet over the course of seven days. 

In this paper, the impact of Shoemaker-Levy 9 on Jupiter is considered as a case study for assessing 
the observability of impacts on giant exoplanets. The comet impact caused visible changes in Jupiter’s 
atmosphere in two major ways on two different timescales. The immediate impacts were very bright as the 
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bolides entered the Jovian atmosphere, but these effects lasted only seconds. Afterward a brightening due 
to thermal emission occurred that lasted on the order of several hours (?). These short duration phenomena 
are unlikely to be observed unless the timing of the comet impact is known in advance, as for SL9, or a 
planet were being continuously monitored. The second major effect of SL9 was to deposit high-altitude 
particulates that blocked sunlight from reaching the deeper atmosphere. The smallest particulates (<0.5 
/rm) stayed at those altitudes for at least two years (?). 

Jupiter’s atmosphere contains 0.2% of methane, which has deep absorption bands in the near-infrared, 
especially near 2.3 /rm in the K-band and 3.3 /im in the L-band (?). When SL9 debris at higher altitudes 
reflected photons from the sun and prevented them from being absorbed by the methane at lower altitudes, 
Jupiter appeared brighter near the impact spots long after the initial effects of the thermal heating had 
ended (?). Because of the tidal disruption, the fragments hit at a wide range of Jovian longitudes, efficiently 
covering the atmosphere at ^45° South latitude. 

For comparable effects to occur with impacts of comets and exoplanets, those planets must have a 
source of opacity in their atmosphere, like methane, whose absorption could be blocked by high-altitude 
particulates. The results from this paper should be applicable to exoplanets as current data and models 
indicate that exoplanets will have such a source of opacity (?). 

2. Observations 

Near-infrared observations of Jupiter were obtained at the Hale 200-inch Telescope at Palomar Observa¬ 
tory on 1994 July 22-27 UT, i.e. the nights following the last of the fragment impacts on 22 July. Followup 
observations were obtained over the next month on 16-19 August and 27-29 August. The only photometric 
nights were July 22-23, July 25-26, and August 28-29 UT. 

Images were taken with the D78 Cassegrain camera, which employed a 256x256 HgCdTe (NICMOS3) 
array with a pixel size of 0.125" square, for a field-of-view of 32x32 arcsec. The camera was rotated to put 
the North-South axis of Jupiter along the detector columns. Over the ~1 month of observations, Jupiter’s 
apparent diameter ranged from approximately 37.7 to 34.2" (radius of 151 to 137 pixels) along the equatorial 
axis and from 35.2 to 32" (radius of 141 to 128 pixels) along the polar axis. More of the planet could be 
seen in the late August images than in the July images, as Jupiter was moving away from Earth during 
that time, going from a distance of about 5.24 AU to 5.78 AU. In addition, to avoid cosmetic defects on the 
array — dead pixels which would have prevented photons from the south pole from being counted — the 
planet was generally placed somewhat higher than center. Therefore, Jupiter was too large to completely 
fit on the detector, and a portion of the top of the planet was cut off in every image. The detector also had 
many bad pixels, but fortunately most of them were in the top half of the array and not near where the 
impact spots appeared. These pixels were replaced with the weighted average flux of surrounding pixels, 
which had minimal impact as the total planetary flux was calculated and not the detailed distribution of 
flux. Through July 25, a second instrument, Spectro-Cam 10, was mounted simultaneously with D78, and 
its pick-off mirror vignetted on the top of the D78 images. On July 26, condensation from its dewar dropped 
water onto D78, causing bright spots on the images. Again, while this contaminated the quality of the top 
part of the image, the lower part, with the impact spots, was unaffected. 

The observations are summarized in Table [TJ Sequences of images of Jupiter were collected across the 
K and L-bands. A typical sequence started with an image in a methane filter centered at 2.276 /rm with a 
bandwidth of 0.17 /rm. The sequence continued with images taken in 1% bandwidth continuously variable 
filter (CVF) at eight central wavelengths in the K-band (2.00, 2.05, 2.10, 2.14, 2.17, 2.20, 2.25, and 2.35 
/im). At the end of the K-band CVF sequence, another methane filter image was taken. The sequences 
finished with L-band CVF observations at 11 central wavelengths (3.00-4.00 /im at 0.10 /im intervals). In 
both bands, individual frame integration times were either 1 or 0.4 seconds, repeated and co-added to reach 
at least 10 seconds total integration time at each filter. In each sequence, the telescope was nodded to sky 
for an image after every image of Jupiter and these sky images were used for background subtraction. 

Generally, two sequences were completed each night. Each image sequence lasted 20 minutes. During 
this time Jupiter rotated ^12 deg. The 2 /im data were previously used in ? for a study of the aerosol 
clouds formed by the impacts and the Jovian winds that dissipated them. 
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date 



times 

UT 

wavelengths 

int time 

sec 

airmass range longitude range 
degrees 

July 

22 

23:14 

to 

23:17 

2.276 

10 

2 

1-2 

July 

22 

23:18 

to 

23:30 

2.0-2.35 

20 

1.9-2 

3-10 

July 

22 

23:37 

to 

23:59 

3.0-4.0 

40 

1.7-1.9 

14-27 

July 

23 

00:00 

to 

00:22 

3.0-4.0 

40 

1.6-1.7 

28-40 

July 

23 

02:14 

to 

02:16 

2.276 

1 

1.4 

114 

July 

23 

02:27 

to 

02:47 

3.0-4.0 

40 

1.4-1.5 

117-130 

July 

23 

02:54 

to 

02:56 

2.276 

20 

1.5 

134 

July 

23 

02:56 

to 

03:04 

2.0-2.35 

20 

1.5 

135-140 

July 

23 

03:06 

to 

03:10 

2.276 

20 

1.5 

141 

July 

25 

23:30 

to 

23:35 

2.276 

10/20 

1.8 

102 

July 

25 

23:35 

to 

23:46 

2.0-2.35 

20 

1.7-1.8 

105-111 

July 

25 

23:47 

to 

23:48 

2.276 

10 

1.7 

112 

July 25- 

■26 

23:49 

to 

00:23 

3.0-4.0 

80 

1.5-1.7 

113-133 

July 

26 

00:41 

to 

01:16 

3.0-4.0 

80 

1.4-1.5 

145-166 

July 

26 

01:30 

to 

01:32 

2.276 

0.8 

1.4 

174 

July 

26 

02:18 

to 

02:19 

2.276 

20 

1.4 

203 

August 

28 

23:02 

to 

23:23 

3.0-4.0 

40 

1.5 

177-187 

August 

28 

23:24 

to 

23:25 

2.276 

20 

1.5 

191 

August 

28 

23:26 

to 

23:36 

2.0-2.35 

20 

1.5 

192-198 

August 

28 

23:37 

to 

23:42 

2.276 

10/20 

1.2-1.5 

199 

August 

29 

00:45 

to 

00:56 

2.276 

10 

1.5 

246 

August 

29 

00:57 

to 

01:07 

2.0-2.35 

20 

1.6 

247-253 

August 

29 

01:08 

to 

01:09 

2.276 

20 

1.6 

254 

August 

29 

01:09 

to 

01:45 

3.0-4.0 

80 

1.6 to 1.7 

255-275 

August 

29 

01:47 

to 

01:48 

2.276 

20 

1.7 

278 

August 

29 

01:49 

to 

01:59 

2.0-2.35 

20 

1.7-1.8 

278-284 

August 

29 

02:00 

to 

02:01 

2.276 

20 

1.8 

285 

August 

29 

02:06 

to 

02:21 

3.0-3.7 

40 

1.9-2.0 

289-297 


Table 1: Sequences of images of Jupiter taken at the given dates, times, and wavelengths. 


star 

K band magnitude 

L band magnitude 

HD 129655 

6.69 

6.66 

HR 5107 

3.06 

3.05 

HR 6707 

3.20 

3.15 

HR 6147 

2.26 

2.18 

HD 161903 

7.02 

6.98 


Table 2: The standard star magnitudes for HD 129655, HR 6707, and HD 161903 were adopted from ? and HR 5107 and HR 
6147 from CIT unpublished measurements (G. Neugebauer, personal communication). 
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Standard stars with known magnitudes (Table ^ were observed on the three photometric nights, as 
summarized in Table |31 These bright stars were nodded on the detector for sky subtraction. HR 5107 was 
streaked with the chopping secondary to avoid saturation. 

date times (UT) stars airmass range 

July 23 03:11 to 03:43 HD 129655 T3“ 

July 26 01:33 to 02:07 HD 129655, HR 5107 1.2-1.3 

August 28-29 23:57 to 00:39 HD 129655, HR 6707, HR 6147, HD 161903 1.2-1.7 

Table 3: Standard stars used on the given nights. 


3. Data Analysis 

3.1. Image Processing 

Basic imaging processing steps were: a linearity correction for the response of each pixel, sky subtraction 
using the offset sky image, and flat-fielding to correct for pixel-to-pixel response variations. The flat-field 
image was created by dark subtracting average sky images at 2.27 /rm and 3.5 ^m for the K-band and 
L-band observations respectively. 

A large number of bad and intermittently bad pixels had to be identified and corrected. To determine 
the location of these pixels, about 50 images with high signal-to-noise ratios were randomly chosen. In each 
image, obviously bad pixels, those with values over two orders of magnitude different from any surrounding 
pixels, were identified first. Then, pixels off by one order of magnitude or more were added to a second list 
in order to cross check. If a pixel was listed for 75% of the files, it was included in a general bad pixel mask. 
To determine which pixels were bad for each specihc image, a program was created that went through each 
file and determined which pixels were not within two orders of magnitude of the median of the surrounding 
pixels. This was combined with the general bad pixel mask to create a mask for each individual Hie. 

Sky subtraction using the offset sky images was imperfect due to the time difference between the image 
of Jupiter and its corresponding blank sky image. An additional sky correction was made by finding an 
average value of the pixels outside the planet by at least 10 pixels. Not included were pixels contaminated 
or vignetted or marked as a bad pixel by the image’s bad pixel mask. After an initial sky average was 
calculated, any pixel whose value was not within 3 sigma of this average was excluded, and a second average 
was calculated. This second average was what was subtracted from the image. 

Each pixel marked as a bad pixel needed to be corrected. Pixels within three pixels of the bad one are 
used to calculate a new value based on their average, with pixels closer to the bad pixel being weighted more 
than pixels further away. Other bad pixels within the three pixel radius are excluded. 

3.2. Photometric Calibration 

The standard star magnitudes for HD 129655, HR 6707, and HD 16190 were adopted from ? and HR 
5107 and HR 6147 from GIT unpublished measurements (G. Neugebauer, personal communication) (see 
Tabled]). Standard stars were observed at airmasses of 1-1.8. 

Basic image processing for the standard star frames was the same as described in the previous section. 
In each of the two standard star images, a 2D Gaussian fit was used to determine the center of the star. 
Aperture photometry was done in a 20 pixel radius except for the elongated images of HR 5107, in which a 
rectangular mask was used. The aperture was large enough to encompass all of the stellar flux. 

The airmass correction was calculated using the standard airmass correction for Palomar. 

3.3. Whole Planet Photometry 

In order to compute the total brightness of Jupiter at each wavelength and time, we had to identify the 
pixels belonging to the planet in each image. Because Jupiter bulges at the equator, its shape was treated 
as an ellipse with values of its semi-major and semi-minor axes required. The equatorial radius of Jupiter 
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Figure 1: An example of a recreation of Jupiter, using an image taken on July 25 at 2.27 pm . The blue box indicates what was 
in the original image. Everything outside of the blue box had to be extrapolated. The sides that were cut off were extrapolated 
by averaging the flux for that latitude. The north pole was extrapolated by duplicating the south pole but at half its flux. The 
scaling in this image is the square root of the flux, so that the different latitudes can be more easily seen. 


is 71492±4 km, and the polar radius is 66854±10 km (?). As the uncertainties in the radii are negligible, 
they were ignored in the calculations. As observed, the radius in pixels was: 

206265 

where R is the radius of Jupiter in km, d is the distance from Earth to Jupiter in km, 206265 is the number 
of arcseconds in one radian, and 0.125 is the pixel scale of the camera. 

An ellipse with those measurements was manually aligned with each image, and the center of the ellipse 
taken as the center of the planet. 

The JPL HORIZONS (?f] ephemeris was used to determine the central longitude of the planet and its 
radius at the time of each image. The planet was divided into three regions based on latitude. The first was 
the region where the planet was impacted, which ranged from approximately 60° south latitude to 30° south 
latitude, the south polar region which is everything south of the spots region (i.e. 30° S), and a region near 
the equator that was relatively symmetric with respect to longitude from approximately 15° south latitude 
to the equator. 

Unobserved areas of the planet were extrapolated as follows. The average pixel value in every detector 
row (latitude) sampling the planet was used to extrapolate the image out to the ephemeris radius. For the 
north polar region, the south pole region was duplicated but at half its flux. The flux was halved because 
the north polar region has less total brightness than the south polar region (?). In each sequence of images, 
the flux from each region was summed. 

We did not have observations of Jupiter in the same filter sequence prior to the impacts. However, we 
can construct a model of a pre-impact Jupiter-like planet, i.e. a gas giant planet with much of its flux 
coming from the poles, by combining two polar regions (i.e. duplicating the measured south polar region in 
the north) and a central region, which is everything between 60° north latitude and 60° south latitude. The 
pixels in the central region were set at the average flux from around Jupiter’s equatorial region between 0° 


^http: //ssd.jpl.nasa.gov/?horizons 


5 









Post-impact Pre-impact 



Figure 2: An example of the model giant planets created at 2.27 pm . These model planets were based on Jupiter with much 
of the flux coming from the poles. They are identical in the two polar regions and the central region from 60° north latitude to 
30° south latitude. The pixels in the central region are set at the average flux from around Jupiter’s equatorial region between 
0° and 10° north latitude. From 30° south latitude to 60° south latitude, the flgure on the left has the post-impact planet 
with impact spots, while the flgure on the right, representing the pre-impact planet, has an extended central region. 


and 10° north latitude. A Jupiter-like planet after the SL9 impact consists of the same two polar regions, 
a similar central region from 60° north latitude to 30° south latitude, and the impact spots from 30° south 
latitude to 60° south latitude. 

For each sequence of images, we measure the surface brightness of the planet in each of the three regions 
and construct a simulated planet. The measured pixel values are converted to magnitude by taking their 
logarithm multiplied by -2.5. A specific airmass correction for each region or planet was added to this 
magnitude. 


AMC = z — {am — 1.0) * amc 

where z is the zeropoint for the wavelength, am is the airmass for the image from which the planet or region 
was taken, and amc is the general airmass correction for the wavelength. 

The apparent planet albedo, i.e. the planetary radiance divided by the incoming Solar flux, is given by: 

]_ _ 1 Q0.4(JC-m) 

F n 

where r is the Jupiter-Sun distance in AU, is the pixel solid angle, K is the apparent magnitude of the 
Sun (e.g. ?). 


4. Results 

Exoplanets will not be spatially resolved, so total integrated planet brightness would be the only measur¬ 
able quantity. Brightness is a function of the planet’s distance from the star — which can be estimated and 
corrected for — and the planet’s total albedo. We examined how this total albedo would compare at three 
different times: prior to the impacts, the early days following the impacts, and a month after the impacts. 

To compare what a planet with impact spots would look like to one without impact spots, the ratio of 
I/F from the model giant planets, such as the ones in Figure [U was taken. The impacts would be most 
easily detected at wavelengths where this ratio is very high. The results are displayed in Figure [3] At most 
wavelengths, the difference between the two was negligible. Wavelengths at which the planet is already very 
bright are those without major sources of atmospheric opacity. For example, in the K-band, the increase 
pre to post impacts in brightness is much less obvious at 2.0-2.10 /im than it is at 2.25-2.35 /rm, as shown in 
FigureHl In the region of deep methane absorption centered at 2.3 /rm, the planet became ~2 times brighter 
after the impacts than it was before. However, in the L-band there was no difference, even at wavelengths 
with deep methane absorption (See section 5.1). 
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Figure 3: (Left) The ratio of I/F between the model planets with impact spots to the model planets without impact spots as 
a function of K-band wavelength. At most wavelengths, the difference between the two was negligible, resulting in a ratio near 
1. However, near 2.3 fira, where there is deep methane absorption, the change was much more noticeable. (Right) Same ratios 
for the L-band. Despite a deep methane absorption band near 3.3 /rm, the ratio shows no change as a function of wavelength. 
Thermal emission leaking out of the planet apparently swamps the signal of the impact clouds. The data from the different 
days are artificially off-set in wavelength for clarity. 
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Figure 4: Reproductions of Jupiter from August 29 at 2.0, 2.1, 2.17, and 2.25 pm. At shorter wavelengths, the impact spots 
are more difficult to see, because the rest of the planet is so bright. In the methane band, which is prominent at the longer 
wavelengths, the rest of the planet is fainter, so the spots are clearer. 
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Figure 5: IjF spectra at K and L bands compared across three dates from July - August 1994. There was very little overall 
change in brightness from July 23, the day after the last of the fragments impacted Jupiter, to August 29, the last day of our 
follow-up observations. Over time, winds dispersed the material that was concentrated in the area of the impacts spots in July. 
So in August, while the individual spots are not as bright, they cover a larger area and the total flux is about the same. The 
data from the different days is artificially off-set in wavelength for clarity. 
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Figure 6: Jupiter in late July versus late August. Over time, winds dispersed the material that was concentrated in the area 
of the impacts spots in July. Jupiter appears smaller in the August figure as it was further away from Earth at the time. The 
scaling in this image is the square root of the flux, so that the different latitudes can be more easily seen. 
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Figure 7: Longitude dependence of J/F at 2.276 pm and at 3.4 pm. The difference in flux is due to the fact that the impacts, as 
well as other atmospheric features such as clouds and storms, were not homogeneously distributed across different longitudes. 


Also of interest was how the total flux changed over the first month following the impacts. Figure [S] 
shows the total brightness of the impacted planet with wavelength in the K-band and L-band compared 
across the photometric days. There was very little overall change in brightness from July 23, the day after 
the last of the fragments impacted Jupiter, to August 29, the last day of our follow-up observations. This 
was also noted in ?. Over time, winds dispersed the material that was concentrated in the area of the 
impacts spots in July. So in August, while the individual spots are not as bright, they cover a larger area 
and the total flux is about the same, as shown in Figure [51 Evidently the settling time of the particulates 
exceeded a month. Thus an observer of Jupiter as an exoplanet would see a brightening of the planet over 
the week-long timescale of the impacts followed by months of constant to very slowly decreasing brightness. 

One complication in deriving a global behavior from our data is that each image sequence covers a 
different range of Jovian longitudes. Since the impact spots were not homogeneously distributed on the 
planet, images centered at different longitudes have different flux as shown in Figure [71 This is also the 
cause in the difference in I/F within one night for a given wavelength, as Jupiter rotates very quickly, with 
a rotation period of just under 10 hours (?, and references therein). To examine the variability of the planet 
with time after the impacts, we use only data taken at longitudes between 165° and 225° to get a more 
direct comparison. These longitudes were chosen because there were data from both months in the K-band. 
Figure [S] shows the ratio of the model planets with impact spots to those without impact spots from only 
those longitudes, and the effect is more obvious. However, the observations in the L-band had no overlap in 
longitude for the photometric days, as shown by Figure [71 This made it impossible to do a direct comparison 
between days. 

5. Discussion and Conclusions 

5.1. Analysis and Explanation of Results 

At the wavelengths where methane has strong absorption features, Jupiter’s atmosphere has very high 
opacity. Results show that the typical atmospheric pressure level probed by K-band observations drops 
from a few hundred mbar outside the 2.3 fim methane band to ^40 mbar inside it (?). The SL9 impact 
spots, which are probably composed of sub-micron aerosols (?), are bright in the near-infrared because they 
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Figure 8: The ratio of I/F between the model planets with impact spots to the model planets without impact spots for images 
with central longitudes between 165° and 225°, comparing I/F from July and August. The data from the different days is 
artificially ofT-set in wavelength for clarity. 


efficiently reflect sunlight before it reaches altitudes where it can be absorbed by methane. Their near- 
infrared brightness is comparable to that generated by the polar hazes, which are composed of hydrocarbon 
aerosols, i.e. submicron particles, generated in Jupiter’s atmosphere by energetic particles launched off lo 
and carried to Jupiter along its magnetic Held lines (??). In general, it is not known if exoplanets will 
have bright poles due to high altitude hazes, such that viewing angle will also affect the contrast between 
comet-deposited material and intrinsic haze. After a month, particles smaller than 0.1 /im would not have 
settled below 40 mbar while particles smaller than 0.5 /im would not settle below 100 mbar (?). 

Although Galileo data shows that the 3.3 /im methane band also probes only the very uppermost layers 
of the atmosphere ~1 mbar, (?) the L-band images of Jupiter show that there are bright bands of emission 
that do not correspond to the reflective areas at K-band. Figure IHl shows no difference between the model 
planet images at that wavelength, so we conclude that there is thermal emission leaking out, comparable in 
magnitude to the scattered light off the impact spots. 
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Figure 9: Reproductions of Jupiter from August 29 at 3.2, 3.4, 3.7, and 3.9 fim. At these wave lengths, the central latitudes 
are much brighter, thus making the spots undetectable, even in the methane band. 
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5.2. Imaging Exoplanets 

At present, the few directly imaged exoplanets known (e.g. HR 8799 bcde (?)) are young, warm bodies 
glowing in their own thermal emission. Although these planets are embedded in massive debris disks, their 
internal emission would likely wash out the effects of impacting planetesimals. Current extreme adaptive 
optics systems designed for planet-finding, such as GPI (?) and SPHERE (?) should directly image young 
luminous planets. The ability to detect cool planets, such as Jupiter with a surface temperature of ~130° K 
(?), is probably left to future giant telescopes and high contrast space missions. Proposed coronagraphs that 
can achieve a contrast ratio of better than 10“® will be able to image these planets with separations of greater 
than 0.2"(?). Atmospheric properties of warm planets are determined by measuring their spectra at multiple 
wavelengths, and the presence of methane and other molecules are used to determine their temperatures and 
compositions. This technique will not change in the era when cooler planets can be observed in scattered 
light. Already, there are ~10 planets detected with the radial velocity method around main sequence stars 
with ages comparable to that of the Sun that are also amenable to high contrast imaging observations. More 
will be found as RV surveys continue for longer durations, and these will be the set of planets most amenable 
to direct imaging with future instruments. (W. Traub, personal communications) 

5.3. Planetesimal Impact Detection 

Detecting these planetesimal impacts is important for several reasons. The effect of planetesimal impact 
spots would be a source of signal and a source of obscuration. Because the planetesimal does not affect 
every longitude equally, bright spots that last many months could allow the rotation rate of an exoplanet 
to be determined using the photometric variations as a function of longitude. However, as accurate spectra 
are necessary to learn important information about planets, anything that might alter the spectra has to 
be taken into account, and given the length of time (i.e. many months) it takes the spots to fade, reflected 
light from the impact spots could obscure the molecular absorption bands normally used to determine 
atmospheric composition. The equivalent widths of molecular absorption lines will be lower than expected 
if high altitude clouds conceal the main sources of atmospheric opacity. Additionally, high impact rates 
could affect the atmospheric composition of giant planets over time, both due to the energy released by the 
impactors and by introducing, and keeping in non-equilibrium, molecular species that would not be found in 
equilibrium (?). For example, Jupiter has water in its stratosphere, i.e. above its water condensation level, 
only because of the SL9 impact (??). While detailed spectroscopy of extrasolar planets such as that of solar 
system planets will never be feasible, knowing the cometary impact rate would allow realistic models of the 
evolution of an exoplanet’s atmosphere. 

To show the same effect that SL9 caused on Jupiter, an exoplanet would need to show deep methane 
absorption. Methane has already been found in some hot Jupiters (e.g. ?). And fortunately, it seems likely 
that cooler exoplanets will also have methane. All of the outer Solar System planets have higher methane 
mixing ratios than Jupiter (????). As these are the coldest planets imaged, with temperatures down to 
^65° K for Neptune (?), methane could be the primary source of opacity for any exoplanet cooler than 
Jupiter. And models indicate that there is significant methane absorption in the K-band in planets with 
effective temperatures up to at least 800° K (?). Even if there are cold exoplanets where methane is not 
a source of opacity, they are likely to have other molecular species, such as carbon dioxide or ammonia; 
comparable measurements could be made at wavelengths where those molecules have absorption bands. 

A planet covered in haze due to abundant aerosols would have less detectable comet impacts because it has 
less detectable atmospheric absorption lines. Aerosol formation is not well understood and is certainly related 
to poorly understood vertical mixing that either does or does not bring methane above the condensation 
pressure. Essentially what the impact does is mimic the reflective effect of such aerosols. When the depths 
of cool exoplanet molecular lines are eventually measured, comet impacts will be a source of particulates 
that could mimic large-coverage hazes. Looking at atmospheres over several years, to see if the particulates 
settle out or not, and determining if the planets sit in systems with debris disks can help distinguish the 
source of the high-altitude particulates. 

For a planet that was impacted as rarely as Jupiter is now, there is a low probability of observing the 
planet during the time when the impact spots were still visible. Planetesimals as small as 200 m could leave 
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long-standing atmospheric debris, which would make detections possible, and those impact more frequently 
— once or twice per decade in the case of Jupiter (?). But a planet being impacted very frequently would 
also be problematic, as a planet that was impacted at the rate per month would show weaker molecular 
absorption bands than one impacted rarely. 

How often a planet is impacted by planetesimals is intimately related to the architecture of the planetary 
system because the small body populations and planets evolve together. Generally speaking, planets interior 
to massive planetesimal belts, planets close to their stars, and planets with larger gravitational cross sections 
are more likely to be impacted. 

Other planetary systems could have significantly different impact rates. All other things being equal, 
the number of comet impacts should be proportional to the number of comets (?). Planetesimals in other 
systems cannot be detected directly, but the debris dust generated from their collisions can be. The number 
of small dust grains in a debris disk is measured based on the IR excess in the system (e.g. ?). Because the 
lifetime of dust is relatively short, largely due to stellar radiation, the presence of dust implies that there 
must be planetesimal collisions creating that dust. To go from the small dust to an estimate of the number 
of planetesimals is model dependent because the intrinsic size distribution of bodies depends on a number 
of factors including the bodies’ size and strength (e.g. ??). However, systems with more dust must have 
more planetesimals. 

Cold debris disks appear to be approximately equally common for planet-hosting and non-planet hosting 
stars (??), although the known planets are all closer to their star than Jupiter is to the Sun. Detectable disks 
are typically 10-100 x dustier than our Edgeworth-Kuiper Belt (??), with correspondingly larger parent 
body populations. We would expect planets in these systems to be subject to frequent impacts. 

Debris disks decay over time due to collisions that erode the parent planetesimals and radiation effects 
that remove the small dust grains. Therefore, impact rates are likely higher in younger systems. However 
? shows that even 10% of star of age similar to the Sun have disks substantially more massive than our 
Edgeworth-Kuiper Belt. 

Solar system architecture also has a large effect on the number of comets. In our Solar System, the 
comet supply has been severely depleted by giant planet migration (?). But that may not be the case in 
all other systems. The ratio of comets to asteroids among planetesimals is a function of how far the stable 
planetesimal belts are from their stars. In our Solar System, asteroids are from the closer and warmer 
Asteroid Belt, while comets come mainly from the Edgeworth-Kuiper Belt and the Oort Cloud, because 
those are at temperatures at which volatiles - a defining component of comets- freeze (?). The disks that 
we know of around other stars are largely cold and analogous to the Edgeworth-Kuiper Belt, and therefore 
capable of delivering impacts to cold giant planets. Oort Cloud equivalents are not yet detectable. A more 
complex factor is the planetary system as a whole. Other planets can increase or decrease the number 
of cometary impacts depending on their sizes and locations, but without detailed information on specific 
systems, it’s not possible to estimate the effects of different architecture on the impact rate. 

There are many other things that could complicate the detection of a planetesimal impact through this 
method. Eor a planet much hotter than Jupiter, the thermal emission would dominate the K-band, like it 
does the L-band with Jupiter, thus making the effect of a planetesimal impact much harder to detect. A 
large asteroid would likely result in similar results as a comet impact. An impact on Jupiter in July 2009 
was by what is believed to be an asteroid under 1 Km in diameter, yet the effects on the planet were thought 
to be relatively similar to that of SL9, especially at first (???). And while rotation rate could be detected, 
an unusual weather system, like the Great Red Spot on Jupiter, would have a similar effect. 
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